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ABSTRACT 


“ 


During the consolidation phase, reinforcement particles of Metal Matrix 
Composites (MMC’s) tend to be non uniformly distributed. The result is that the material 
properties of the composite materials are not as good as those originally desired. Through 
large amounts of straining, homogeneity can be achieved. Finite element models of 
MMC’s undergoing different thermomechanical processes (TMP’s) to true strains of 
approximately 1.2 were generated. The models consist of particle clusters within the 
particle-depleted matrix. The particle clusters were modeled by either a smeared model in 
which the particles refine the grains in the cluster, or a discrete model of the particles 
within clusters. The smeared and discrete models qualitatively agreed with each other. 
The results suggest that the best TMP to reach a state of reinforcement particle 


homogeneity was a hot worked, low strain rate TMP. 
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I. INTRODUCTION 

Composite materials have been used for centuries. Straw in bricks and bulls hair in 
plaster were used to increase the strength and durability of ancient structures. In more 
recent years reinforced concrete was one example of a composite material which was a 
major breakthrough in building technology (Caddell, R. M., 1980; and Parton , V. Z., and 
Kudryautsev, B. A., 1993). 

A modern definition of a composite would be “the combination of a reinforcement 
material (such as a particle or fiber) in a matrix or binder material” (Strong, A.B., 1989). 
A more refined definition of a composite includes: 1) The material must be man made. 
2) It must be a combination of at least two chemically distinct materials with a distinct 
interface separating the constituents. 3) The constituent materials forming the composite 
must be combined three dimensionally. 4) It should be created to obtain properties that 
would not otherwise be achieved by any of the individual constituents (Taya, M., and 
Arsenault, R. J., 1989). Under this definition the first modern application of a composite 
material was a glass reinforced phenolic-nylon fishing pole in 1945. The first production 
application of Metal Matrix Composites (MMC’s) was boron/aluminum tubes in 1974 as 
components of the space shuttle. Over the past few decades the manufacture and use of 
composites have exploded as engineers and consumers have kept looking for stronger, 
lighter, and more compatible materials. Examples of the application of these improved 
composite materials include automobile engines, turbine engine components, and 


aerospace structural components in which the material must operate effectively in harsh 


environments, and higher strength to weight ratios are critical. The research into new 
composite materials as well as manufacturing techniques has been of major interest and 
will continue to be so in the future as engineers strive to design =e and more durable 
systems. 

Particle-reinforced metal matrix composites (MMC’s) consist of relatively small, 
hard particles embedded within a softer metal matrix. The reinforcement particle volume 
fraction may range from 5% to 40 %. Particle sizes are typically one micron or greater. 
The nature and function of the particles can vary widely, resulting in significantly different 
material properties. Particulate composites differ from flake and fiber reinforced types, 
and the additive constituent 1s usually random tn distribution. Because of this, particulate 
composites are usually assumed to be isotropic (Schwartz, M. M., 1992). Although 
particle reinforced MMC’s exhibit higher strength, stiffness and wear resistance, their use 
has been limited due to relatively low ductility and fracture toughness. 

In order to achieve the improved mechanical properties and an acceptable fracture 
toughness desired in a particulate MMC, the particles must be dispersed uniformly and 
homogeneously throughout the matrix material. The production of a MMC consists of 
two major steps: the initial consolidation of the material, followed by the forming of the 
material into a useful shape. The fabrication, using either casting methods or powder 
metallurgy techniques, generally does not result in a uniform dispersal of the particles in 
the matrix (Taya, M., and Arsenault, R. J., 1989). Therefore, regardless of the 
consolidation/fabrication technique, the reinforcing particles of MMC’s tend to be non- 


uniformly distributed. As a result, the material properties of such composite materials are 


poor. Because of this the materials may be subjected to large amounts of straining in the 
forming process in an effort to achieve a uniform and homogeneous microstructure. 
McNelley and Ballou (1995) performed an experimental feds of the mechanisms 
by which this particle redistribution was taking place in a 6061 AL-AL,0; MMC. Their 
study examined the role of temperature during the deformation process and the 
experimental results indicated that the particle distribution became more uniform using a 
hot thermomechanical process (TMP) as compared to using a cold TMP. 

The purpose of this study was to model, using the Finite Element Method (FEM), 
the hot and cold TMP’s of particle-reinforced MMC’s. From the FEM models, a 
qualitative analysis was performed to determine the effects of temperature as well as strain 
rate in determining the ability of a cluster of particles, or the particles themselves, to move 
within the matrix and thus achieve a state of homogeneity while undergoing large ene 
A qualitative evaluation was performed of the stresses and strains achieved within the 
matrix while undergoing the TMP’s. The objective was to then determine, from a 
qualitative viewpoint, what condition of TMP’s would result in, as close as possible, the 


homogeneous particle redistribution to achieve the desired mechanical properties of a 


particle-reinforced MMC. 





Ir. BACKGROUND 


A. THEORY OF CONSTITUTIVE EQUATIONS * 

The theory of plasticity 1s concerned with the study of stresses and strains in 
plastically deformed materials such as metals. The theory is used to develop constitutive 
laws relating stresses to strains in plastically deformed solids (Hill, R., 1956). These 
constitutive laws may then be used to model large deformations of the solid material while 
it 1s either subjected to cold-working or hot-working. In the present discussion of the 
constitutive laws only isotropic material 1s considered. 

1. Strain Hardening Plasticity 

Plasticity of a cold material occurs when the applied stress acting on the material 
Betis it yleld strength (1.e. 0 > Oo; where o = applied stress and o 9 = yield strength). 
Once this occurs, plastic flow begins. Below in Figure 1 is a plot of stress vs. strain 
depicting the behavior of a material that is plastically deformed. The solid line from A to 
C shows the loading of the material to an applied stress above yield strength with 
subsequent unloading along the dashed line from C to D. The subsequent behavior of the 
material following this process is such that the yield strength will increase to the value of 
the maximum applied stress during this process, and the elastic modulus of the material 
will remain constant. This phenomena is known as strain hardening. 

There are two commonly used general theories which predict initiation of plastic 


flow, Tresca’s flow criterion and von Mises’ flow criterion. 


Stress vs. Strain 





Figure 1: Stress vs. Strain behavior for a material which ts plastically deformed. 
O app = applied stress and o q = yield stress. 
Tresca’s flow criterion is based on the maximum shear stress being greater than the 

shear yield or flow strength of the material. Specifically, 

Tes = IG (235 
where 7 o = Shear yield stress, and the maximum shear stress 1s expressed as 

Grane Cite Gi) 2 22 
Therefore, Eq. (2.1) can be rewritten as 

O 1-0 3 16 6 (2.3) 
where 0 1}2> 02203, and 01,02, ando;3 are the principal stresses and 

To= Oo/2 (2.4) 
Plastic flow occurs when Tmax > To, For a biaxial state of stress (i.e. one of the principal 
stresses are zero, o 3 In this case) the following plot in Figure 2 1s generated. This 1s a two 


dimensional plot of the yield locus. The area inside the yield locus represents the principal 


stress combinations in which no yielding occurs, while stress combinations outside the 


yield locus cause yielding and plastic flow begins. In three dimensions this plot look like a 


™~ 


hexagonal prism with its central axis aligned 45° from each of the principal axes. 


Tresca's Yield Locus 





Eicure 2: Tresca’s Yield Locus. 


The von Mises’ flow criterion 1s based on distortion energy. Yielding begins when 
the elastic energy of distortion reaches a critical value. Hydrostatic pressure is assumed 
not to contribute to yielding because it only produces a dilatational strain energy (Hill, R., 
1956). Asa result, yielding begins where o > oO 9 where 
o et = (I/V2 )[(o1- 02) + (02-03) + (3-01) }”. (2.5) 
Again, for a biaxial stress state it can be shown that 


a el es eae) = aes (2.6) 


and the yield locus for the von Mises’ criterion 1s shown below in Figure 3. For three- 
dimensional principal stresses, the shape takes the form of a right circular cylinder with its 
central axis aligned 45° from all of the principal axes. 


It should be noted here that for isotropic materials and isotropic strain hardening, 


the nght circular cylinder will expand equally in three dimensions as o 9 increases. 


von Mises Yield Locus 





Figure 3: The von Mises Yield Locus. 


In terms of the stress tensor, o ; , the following may be shown. First 
Gir Oat Oo (Qa 
where o’ jj 1s the deviatoric stress, 64; is the Kroneker delta and o , is the hydrostatic 
stress component. Plastic flow is only dependent of the deviatoric stress component. For 
a state of uniaxial stress, equations 2.8 and 2.9 may be derived: 


G 1--26 >, (em (23) 


O€ } = .-2 O€ 2 = -2 OE€ 3 (2.9) 


Equation (2.9) assumes that no volume change occurs during the plastic deformation. 


e 


Therefore the following equations hold true: 


O€ \/ 0&2 = -2 (2.10) 
661562 = 0’ lo"? (2.11) 
OEY 1 — O€2/0 2 = 0&3 0" 3 (2.12) 


Based on the preceding equations, it can be assumed that 

be i/o' i= OA (2.13) 
where 6A is the scale factor of proportionality. This results in the Levy-Mises equation 
which is rewritten as 

6&5, = o' OA. (2.14) 
The extension of the Levy-Mises equation to account for the elastic component of strain 
was performed by Prandtl and Reuss and takes the form, 

66°; = o'3(3 60 (20 H’) (2215) 
whereo is the effective von Mises stress and H’ is the slope of the equivalent 
stress/plastic strain curve (Hill, R., 1956). However, 

(OoGiH')=<6 6,5 (2.16) 
where 6 €, 1s the equivalent plastic strain increment which in turn is defined as 

5 Ep = ((2/3) de"; Se)" (2.17) 
Thus the Reuss-Prandtl equation becomes 


de; =o ACs oe p OE » (2.18) 


It is the Reuss-Prandtl equation which relates the incremental plastic strain to the 

deviatoric stress component, the instantaneous hardening rate and the effective von Mises 
~ 

stress. Equation (2.18) is the constitutive equation relating stress and strain which is used 

for the cold worked models in this work. 

Dn Creep 

The behavior of the solid at elevated temperatures 1s governed by thermally 
activated processes. The strength of solids at high temperatures is strongly dependent on 
strain rate, €, and temperature. Strain hardening is reduced because recovery processes 
predominate and the strength of the material can be characterized according to its creep 
resistance. The two environmental factors influencing creep are temperature and stress 
(Sherby, O. D., and Burke, P. M., 1968). 

Figure 4 shows the three stages of creep for a polycrystalline metallic solid. The 
first stage is known as primary creep; in this interval of strain the metal is hardening, and 
this 1s apparent in the decreasing creep rate. The second stage of creep, also known as 
secondary creep, 1s the stage in which the creep rate is constant. In this stage the 
substructure developed has a strong dependency on stress and temperature. The third 
stage, known as tertiary creep, is often associated with localized necking and void 
formation under tension. This stage leads to material failure, and is characterized by an 
increasing strain rate (Sherby, O. D., and Burke, P. M., 1968). As either temperature or 
stress increases, the slope and the strain rate of all three regions increases. The secondary 


stage of creep where the creep rate in constant is the area of interest for hot working. 


10 


Stages of Creep 





Time 


Figure 4: The three stages of creep shown as Strain, €, vs. time. 


Creep is a thermally activated process. Above temperatures of 0.5T,,, where Tm 1S 
the melting temperature of the material, the creep rate 1s dependent on dislocation climb. 
This is often referred to as dislocation creep. Because climb ts a diffusion-controlled 


process the steady state creep rate tends to be proportional to the diffusion coefficient: 


g€ =kDf(o ) (2.19) 
where 
D = Doexp(-Q/RT) (2.20) 


is the diffusion coefficient, T = absolute temperature, R = universal gas constant, 
Q = activation energy for creep, and Do and k are constants. This indicates a direct 


relationship between € and D. Because of this direct relationship it can be shown 


E = K;f(o) (2221) 
where K;=kDo 
or 


€ =K;o". (222) 


At low stresses creep 1s believed not to be due to dislocation climb but due to stress 
“e 

directed atom migration. Here, the stress will bias the direction of vacancy diffusion 

within grains and they will elongate in the direction of the applied tensile stress; this is 


often referred to as diffusional creep or as the Nabarro-Herring creep mechanism as 


illustrated below in Figure 5. 


he, sgl 
- Cee c 


Ficuke: Schematic of the Nabarro-Herring diffusional creep model. 
“The vacancy concentration in region | is higher than the 
vacancy concentration tn region 2 for the central grain illustrated. 
Vacancies will therefore flow in the direction shown by the arrows, or, 
which 1s the same, atoms will flow in the opposite direction leading to 
elongation of the grain in the tensile direction of the applied stress” 
(Sherby, O. D., and Burke, P. M., 1968). 


The type of creep relation given by Nabarro-Herring can be written as 

on (2,239) 
and is thus known as the linear law. At higher stresses the intermediate stress law 1s where 
dislocation climb or the non-conservative motion or jogs of screw dislocations 1s the 


mechanism of creep and expressed as 


é=Kjo° (2.24) 
for pure metals (Sherby, O. D., and Burke, P. M., 1968). 
If creep 1s dependent on two thermally activated processes, 1.e., dislocation creep 


and diffusional creep, each contributing independently to plastic flow, the fastest process is 


always rate controlling. This 1s shown below in Figure 6. 





a Process 2 
\ 
\ 
\ 
= Constant , Process 1 
ape ay ree ee Na 
Wicker 
Figure 6: Schematic of two different processes of creep. “When two independent 


processes control creep, the fastest 1s rate controlling. If each process has a 
different activation energy the observed activation energy for creep will 
change with temperature, increasing with increasing temperature” (Sherby, 
O. D., Burke, and P. M., 1968). 
Grain Boundary Sliding (GBS) is also a form of diffusional creep where the grains slide 
along the boundaries after diffusion of vacancies within the grains has occurred. 
Normally in the intermediate stress range it 1s not a dominant mechanism due to the 


relatively large grain sizes found in pure metal. GBS typically becomes the dominant 


mechanism of creep for grain sizes less than ten microns. For GBS the strain rate-stress 


relationship has been experimentally shown to be (Sherby, O. D., and Wadsworth, J., 
1982) 
é=K307. | (2:25) 

Because of the high density of particles within the cluster regions of MMC’s it is 
thought that Particle Stimulated Nucleation (PSN) of recrystallization within clusters of 
reinforcement particles will result in different elevated-temperature flow characteristics 
within these regions. These regions exhibit recrystallization in the early stages (0.1-0.25 
true strain) of processing and thus behave as fine grained structures. Fine-grain metals 
may deform by grain boundary sliding (GBS) at elevated temperatures and are weaker 
than coarse grained metals. If GBS is the dominant deformation mechanism within 
clusters than they may deform at different rates when compared to the matrix. Such 
differences in flow behavior are analyzed _ 

The constitutive laws thus developed for hot worked models are shown in 
equation (2.24), € =K2o°, for the matrix material; and in equation (2.25), é = K30°, 
for the material within clusters if recrystallization has occurred. 

B. FINITE ELEMENT METHOD (FEM) 

The finite element code used to perform this work was LS-DYNA3D, a vectorized 
explicit three-dimensional finite element code for analyzing the large deformation dynamic 
response of inelastic solids. The code was developed by the Livermore Software 
Technology Corporation (LSTC) in Livermore, California (Hallquist, J. O., 1991). The 
FEM program is based on the Lagragian description of the body. The central difference 


technique is also utilized in the code for time integration. For further technical and 


14 


theoretical information about the program the reader is referred to the theoretical manual 
(Hallquist, J. O., 1991). In the models eight-node hexahedron solid elements were used 
under plane strain conditions. The detailed description of the pee, including the 
meshes and boundary conditions as well as the constitutive equations, is provided in 
subsequent chapters. 

In the present model each cluster zone remains as a single entity even after large 
deformation. In other words, the model does not allow for breaking down of a single 
cluster into smaller clusters and for mixing within the surrounding matrix material. This 
mixing phenomenon certainly happens in actual deformation processing of composites and 
further experimental study is needed to gain better insight into the actual physical 
mechanisms occurring within MMC’s during TMP. The subsequent study should improve 


the modeling technique to more closely describe the actual phenomena. For this purpose a 


Eulerian-Lagrangian technique of the FEM will likely be more suitable. 


le 





Il. SMEARED PARTICLE CLUSTER MODELS OF METAL MATRIX 
COMPOSITES 


otk 


A. DEVELOPMENT OF THE MODELS FOR COLD WORKING 

A few assumptions were made in the development of the particle reinforced MMC 
models. In the initial consolidation of a particle reinforced MMC the reinforcement 
particles are assumed to be non-homogeneously distributed within the matrix material, 
there being local particle-rich areas interspersed with matrix-rich areas. This circumstance 
was simplified in the model; two separate homogeneous entities were assumed: one is the 
particle-rich area, modeled as a cluster, and the other is the matrix-rich area modeled as 
pure matrix material. The shape of the cluster regions was modeled as a cylindrical 
region, or a circular area in two-dimensional space. This can be seen in Figures 7 to 10 
that show the different meshes prior to deformation. 

An assumption made specifically for the cold deformation model was that the 
material behaved as an elastic-plastic material with isotropic strain hardening. The related 
theory of plasticity was discussed in Chapter II. Table 1 lists the material properties used 
for the cluster material and matrix material used in the cold rolled models. The cluster 
region is modeled as a homogeneous material with the effects of the increased fraction of 
reinforcement particles within the cluster being increased Young’s modulus, E, and the 
strain hardening coefficient to 150% of their values compared to the pure matrix material. 
Such an increased modulus E, corresponds to a SiC or AL2O3 particle volume fraction of 


0.25 within the clusters residing in an Al matrix. Because no particles are assumed to be 


present in the matrix, the modulus, E, is that of unreinforced Al. A similar proportional 


increase in strain hardening rate was assumed for the purpose of this analysis. 
~. 
Table 1: Material properties of the matrix material and cluster material used in the 
cold worked models. 
Material Properties of the matrix material and cluster material 


Material type Matrix material Cluster Material 


Yield Stress AO ksi 45 ksi 
Young’s Modulus, E 10.4*10° psi 15.6*10° psi 
Strain Hardening Coeff. 100 ksi 150 ksi 





B. DEVELOPMENT OF THE FEM MESHES 

The overall geometry of the MMC model consists of a three dimensional block 
containing zones representing the clusters. There are four basic mesh geometry’s shown 
below in Figures 7 to 10. The boundary conditions imposed on the models are as follow 
The symmetry boundary condition 1s imposed on the left hand side of the model where the 
left hand nodes are constrained to allow motion in the vertical direction only. The model 
is also constrained on the front and the back surfaces to correspond to the plane strain 
condition. The bottom of the model is set on a stationary rigid body with no friction 
between the material and the rigid body. This allows the material to slide along the rigid 
body surface. The top of the model is also against a rigid body with no friction, but the 
top rigid body is in motion in the downward direction at a constant velocity. The nght 
hand side of the model is free to OU however, when in contact with the top or bottom 
rigid body there is no friction between the rigid bodies and the model 1s therefore allowed 


to slide. In Figure 7 the model contains two particle clusters located side by side. Figure 


18 


8 shows two Clusters in a vertical orientation with the same spacing as the previous model. 
In this and subsequent figures of the mesh, the two rigid bodies at the top and bottom of 


the block are not shown. Figure 9 shows six clusters in a square arrangement and Figure 


10 shows these same six clusters with identical spacing but in a diamond arrangement. 
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Figure 7: Mesh of two cluster regions oriented horizontally. 





Figure 8: Mesh of two cluster regions oriented vertically. 
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Picupe.3: Mesh of six cluster regions in a repeating square arrangement. 
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Ficure 10) Mesh of six cluster regions oriented in a repeating diamond arrangement. 


Cc. STRESS AND STRAIN RESULTS OF COLD WORKING 

All models were subjected to a large deformation (corresponding to a natural 
strain of about 1.2). The fringe plots of effective von Misses stress and the fringe plots of 
effective plastic strain for all the cold worked models at a state of 1.20 true strain are 
shown in Appendix A; the exception 1s that the diamond shaped model is at a true strain of 
1.13. The models of cold deformation indicate significantly higher stresses within the 
cluster regions and relatively lower strains as compared to the surrounding matrix 
material. Looking specifically at-the diamond cluster arrangement, for example, in Figure 
11 it can be seen that the stresses in the cluster region range up to 233 ksi while in the 
stresses in the matrix region are as low as 102 ksi. 
This 1s expected behavior of all the models because the clusters have higher yield strength 
and stiffer modulus. Looking at the fringes of effective plastic strain, shown in Figure 12, 
it is noted that the higher strain regions in the range of up to 1.40 effective plastic strain 
are within the matrix material at the center of the diamond shaped cluster arrangement. 
The lowest effective plastic strains are experienced in the cluster regions themselves and 
range from 0.65 to 1.22. 
D. DEVELOPMENT OF HOT WORKED MODELS 

The models of hot deformation used the same meshes as the cold work models and 
were shown previously in Figures 7 to 10. This was done so that a direct qualitative 
comparison could be made between the two different TMP approaches. These models 
invoke the power law of viscoplasticity in the intermediate stress range of creep (Sherby, 


O. D., and Burke, P. M., 1968) and utilize the constitutive equation: 





Or 


DIF ZCOLD fringe levels 


time = 3. 59986E~-62 1.48867a5 * 
fringes af eff. stress (u-m) 


min= 1.@81L9E+@5 in element 1999 
max= 2@.331E+@S in element 2728 
1.511F&+85 ; 


1.6226485 ; 
1.733E+85 » 


1. B44E405 > 


1.9566+@5 | 





2. 8676485 
2.178E+@65 


2.289E+85 % 


: 2.408E+85 
Zz x 


LS-TAURUS S38.4 Mau95 


Fringe plot of effective von Mises stress for the diamond distribution model 


after cold working to a true strain of 1.13. 
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Fringe plot of effective plastic strain for the diamond distribution model 


after cold working to a true strain of 1.13. 


Figure 12: 
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Giake &. (3.1) 
In Eq. (3.1) m=0 and n and k vary depending on the location in the MMC. 

The following additional assumptions are used in the high temperature modeling 
process. First, the temperature during the hot working is assumed to be above 0.5T,, 
where T,, 1s the melting temperature of the material. Next, it is assumed the matrix 
material behaves such that the mechanism for creep is dislocation climb and therefore n=5 
in Eq. (3.1). This assumption was discussed earlier in Chapter II. In the cluster material 
grain boundary sliding (GBS) is assumed to be the mechanism of creep. Because of 
particle stimulated nucleation (PSN) of recrystallization within clusters after very small 
Strains, it 1s assumed that grain boundary sliding (GBS) becomes the predominate 
deformation mechanism and therefore n=2 in Eq. (3.1) (Sherby, O. D., Wadsworth, J., 
1982). Because the constitutive law is ace On strain rates, the models were deformed 
using three different nominal ea rates. These strain rates include a high strain rate of 


-j D ay 
_ and a low strain rate of 0.333 5 ae 


33.3 s’. an intermediate strain rate of 3.33 s 
determine the k values to be used in the Eq. (3.1), plots of strain rate divided by the 
diffusion coefficient, D, vs. creep stress were used (Sherby, O. D., Wadsworth, J., 1982). 
The n=2 power law and the n=5 power law were assumed to intersect at the intermediate 
strain rate. This is illustrated below in Figure 13. Note that in Figure 13, the line with the 
steepest slope represents the dislocation climb creep mechanism while the other line 
represents the GBS mechanism. It was then determined from the graph that, at a nominal 


stress of 5000 psi and a nominal strain rate of 3.33 s’, that k = 3937 for the matrix 


material and k = 2732 for the cluster material. Note from Figure 13 that, if the strain rate 





is high, the n=2 power law results in significantly higher stresses. In reality this does not 
occur because the n=5 power law will dominate (i.e. dislocation climb) the physical 


behavior. Therefore the ‘k’ value for the cluster material equals the ‘k’ value for the 


matrix material in the high strain rate model. 
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Figure 13: Log-Log plot of intermediate creep stress vs. (strain rate/D) where D = 
Diffusion Coefficient. 

It may be argued that the Nabarro and Herring (NH) diffusional creep law, where 
n=1 in Eq. (3.1), is the mechanism by which creep would occur in the fine-grained 
clusters. This is not believed to be the case here; because NH creep is normally associated 
only with deformation near the melting temperature (T > 0.9 T,,,) and in the very lowest 
creep stress range (Sherby, O.D., and Burke, P.M., 1968). However, one model run using 
the square cluster arrangement was generated using the linear law, n=] and k = 1501.5 in 


Eq.(3.1). The resulting cluster response was very similar to the results using the n=2 


power law. Because of this, the low and intermediate strain rate models were all 
generated using the n=2 power law. 
ie STRESS AND STRAIN RESULTS OF HOT WORKING 

The fringe plots of effective von Misses stress and the fringe plots of effective 
plastic strain for all the hot deformation models at a state of approximately 1.20 tme 
strain, with the exception of the diamond cluster arrangement which was deformed again 
to 1.13 true strain, are shown in Appendix B. It is noted that in the high strain-rate 
models no distinction between the cluster regions and the matrix material is evident as 
shown in Figures 14 and 15. This result is expected because the mechanism of creep is 
dislocation climb within both the clusters and the matrix. 

In the intermediate strain rate analysis the higher stresses were observed to be 
within the cluster material. The effective ERIE strains are also higher within the clusters 
(note Figures 16 and 17). In the earlier stages of hot working at the intermediate strain 
rate, the models do exhibit lower stresses within the clusters. Not until a value of 1.0 true 
strain is reached do the higher stresses in the cluster region become apparent. Near the 
end of the deformation process the strain rates are higher, because of the prescribed 
constant velocity of the ngid body, and the mechanism of creep within the cluster regions 
should change from GBS to dislocation climb. The computer code allows for the 
constitutive law to be expressed as a power function, 1.e., €; = Kio’ or é2= Koo’, 
where in fact the form of the constitutive law when one has differing grain sizes should be 


ET=E1+E2. GZ 
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Fringe plot of effective von Mises stress for the two cluster vertical 


Figure 14: 


arranged model after hot working at a high strain rate to a true strain of 
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Fringe plot of effective plastic strain for the two cluster vertical arranged 


model after hot working at a high strain rate to a true strain of 1.20. 
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Fringe plot of effective von Mises stress for the diamond arranged model 
after hot working at an intermediate strain rate to a true strain of 1.13. 
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This additive form of the rate law is not accepted by the current computer codes. It is at 
an intermediate strain-rate where this limitation will result in significant inaccuracies and 
cannot be ignored. This is one of the limitations of the current acs and needs to be 
improved in future study. 

The low strain rate models result definitively in lower stresses as well as higher 
effective plastic strains in the cluster material (note Figures 18 and 19). The one exception 
to this is that in the square model the strain in the cluster does not differ from the matrix 


material (see Figure 20). The comparison of the linear law and the n=2 power law models 


at the low strain rate is illustrated in Figures 20 to 23. 
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Fringe plot of effective von Mises stress for the two cluster horizontal 


Figure 18: 


model after hot working at a low strain rate to a true strain of 1.20. 
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Fringe plot of effective plastic strain for the two cluster horizontal model 


after hot working at a low Strain rate to a true strain of 1.20. 
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Fringe plot of effective von Mises stress for the square arranged model 
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after hot working at a low strain rate to a true strain of 1.20. 
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Fringe plot of effective plastic strain for the square arranged model after 


Figure 21: 


hot working at a low Strain rate to a true strain of 1.20. 
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Fringe plot of effective von Mises stress for the square arranged model 


Figure 22: 


using the linear law assumption after hot working at a low strain rate to a 


true strain of ].20. 
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Fringe plot of effective plastic strain for the square arranged model using 
the linear law assumption after hot working at a low strain rate to a true 
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Figure 23: 





|S COMPARISON OF DEFORMATION DURING HOT WORKING AND 
COLD WORKING 


Prior to discussing the results of the FEM analysis it 1s worthwhile to review 
microstructural analysis (McNelley, T. R. and Ballou, M. A., 1995) illustrating the 
different effects of cold working compared to hot working. In Figure 24 a non-uniform 
particle distribution is evident. Particle clusters and a region free of particles are seen in 
the as-cast 6061 AL-10 vol. pct. AL203; MMC. Figure 25 shows the same material after 
hot working to a true strain of 1.10. At this point a relatively homogeneous distribution of 
the particles throughout the matrix material has been achieved. Figure 26 shows the same 
material after hot working to a true strain of 0.51, followed subsequently by cold 


deformation to attain a total true strain of 1.20. Although the particles are bettér 
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Figure 24: An optical micrograph of an as cast 6061 Al] - 10 vol. pct. Al,Os3 illustrating 
a non uniform particle distribution (McNelley, T.R. and Ballou, M. A., 
1995): 





nigube 25: A montage of the 6061 Al - 10 vol. pct. Al,O; MMC after hot forging to a 
true strain of 1.10 (McNelley, T.R. and Ballou, M. A., 1995). 
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Proure 26: A montage of the 6061 Al - 10 vol. pct. Al,O; MMC after hot forging to a 
true strain of 0.51 followed by addition cold deformation to a total strain 


of 1.20 (McNelley, T.R. and Ballou, M. A., 1995). 








dispersed, compared to the as-cast material, they are still not as homogeneously 
distributed as the particles in Figure 25; instead, elongated stringers of particles are seen, 
representing the residual influence of clustering in the as-cast material. 

From such microstructural data McNelley and Ballou (1995) concluded that 
particles were more effectively redistributed during hot deformation when compared to 
deformation at lower temperatures. It was also suggested that recrystallization, i.e. the 
formation of new, strain-free grains, maybe responsible for this effective particle 
redistribution. This phenomenon is evident in Figure 25, where refined grains are seen 
throughout; within clusters, such grains appear relatively refined in size. 

In this attempt to simulate the effect of such microstructural refinement, the 
method used for determining the differences in particle to redistribute was to examine the 
cluster deformation by evaluating the height-to-width ratio prior to and after deformation, 
the amount of vertical compression of the cluster material in relation to the matrix 
material, and the horizontal displacement of the cluster centers during the deformation 
process. If the cluster is more readily deformed and more easily moved about within the 
matrix then it is likely that the particle distribution can also more readily be homogenized 
for such conditions. 

Prior to deformation, all height-to-width ratios of the cluster regions are 1.0. 
Figures 27 to 30 show the height-to-width ratios for the clusters after deformation. Each 
figure applies to one mesh geometry, and the simulations of different TMP conditions 


differ within each figure. In Figures 27 to 29 the true strain is 1.28, while for Figure 30 
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the true strain is 1.13. Following is a description of the abbreviations used in Figures 27 


to 36. 

CW: The cold working process. 

HW-high: The hot working process under a high strain rate. 

HW-int: The hot working process under an intermediate strain rate. 

HW-low: The hot working process under a low strain rate. 

HWL-int: The hot working process using the linear law in the cluster material and an 
intermediate strain rate. 

HWL-low: — The hot working process using the linear law in the cluster material and a 


low strain rate. 
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Ficitnew a. Height to Width ratios for the two clusters arranged horizontally. 
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Height to Width ratios for the six clusters in a square array. 
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Figure 29 
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Figure 30: Height to Width ratios for the six clusters in a diamond array. 


From Figures 27 to 30 the observation can be made that the ratios are definitely 
smaller for the hot working araneeees when compared to the cold working process, 1.e., 
the clusters deform more readily during hot working. This would suggest that the MMC 
would more readily homogenize during hot working. Also, the ratios are also smaller for 
the low strain rate vs. the high strain rate analysis. This further suggests that the MMC 
would more readily homogenize using a low strain rate hot working process. 

Next, a close examination of the vertical compression of the cluster region in 
relation to the matrix material was performed at the centerline of the block. This may 
provide insight into the relative deformability of the cluster and matrix regions. Two 
cluster models were used for this analysis: The vertical and horizontal arrangement of two 


clusters. The analysis are presented in terms of a percentage of cluster material in the 
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vertical direction as compared to the height of the model. The results follow in Figure 31 
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Once again it can be seen from Figure 31 and 32 that the cluster is more readily deformed 
in the hot models vs. the cold model. Again this suggests that the MMC would more 
readily homogenize during the hot working process. 

Now, the displacement of the cluster centers in the horizontal direction is 
compared in Figures 33 to 35. In Figures 33 and 34 trie strain 1s 1.28, and for Figure 35 


true strain is 1.13. The two cluster vertical arrangement is not shown since no horizontal 


displacement of the cluster center 1s allowed due to the symmetry boundary condition. 





Horizontal Displacement of cluster center 
for two cluster horizontal model 





Displacement 


Figure 33: Horizontal displacement of the cluster center for the two-cluster horizontal 
arrangement at a true strain of 1.28. 
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Horizontal displacement of the cluster centers for the square 
model 
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In Figure 34 data for only four clusters are shown. This is because the two clusters at the 

left of the square model have no displacement in the horizontal direction. Figure 35 
~ 

displays data for five clusters in the diamond array for the same reason. It can be seen in 

Figures 33 to 35 that the horizontal displacement is, once again, greater in the hot worked 

models and this again leads to the conclusion that the MMC would more readily 

homogenize during the hot working process. 

All of the previously displayed results also indicate that the deformation associated 
with clusters increases as the strain rate is lowered during the hot worked process. There 
are significant difference between the high strain rate and the intermediate strain rate 
cases, while the differences between the intermediate and low strain rate models is evident 
but smaller. The conclusion is then made, based on the assumptions of this work, that the 
best TMP to follow to reach a homogeneous particle distribution in the least strain is a 


relatively low strain rate during a hot working process. A summary of all the results is 


illustrated in Figure 36 for the two-cluster horizontal array. 
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Visual view of cluster deformiution for the two cluster horizontal models 





Figure 36: Visual view of the cluster deformation after hot and cold working to a true 
strain of 1.28 in the two cluster horizontal models. 

One area of concern in this study was how the free boundary of the mesh affects 
the horizontal movement of the clusters. To answer this question the six-cluster diamond 
array mesh was modified by adding a large amount of matrix material to the right hand 
side of the mesh, thereby extending the mesh in the horizontal direction. This tsolated the 
clusters from the free boundary effects. After deformation the horizontal displacement 
was compared to that of the diamond array without the mesh extension. It was found that 
the clusters did not displace at all during the first 0.02 of true strain as compared to the 
first 0.01 true strain in the non-extended mesh. However, once the clusters began to 


move, both meshes show the same horizontal displacement of the cluster centers. 
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Therefore, from a qualitative point of view, the free boundary has no effect on the 


horizontal displacement or cluster deformation. 
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IV. DISCRETE PARTICLE MODELS OF METAL MATRIX COMPOSITES 


™ 


A. DEVELOPMENT OF DISCRETE MODELS AND FEM MESHES 

The preceding analyses neglected the particles and assumed the only influence of 
the particles was on grain size, during elevated temperature deformation. This ignores the 
effect of the hard particles on stress states and the modulus of the composite. Therefore, 
individual reinforcement particles were modeled discretely in cluster regions of the mesh 
as shown in Figures 37 to 41. The particles are still in clusters within the matrix and 


furthermore arranged in the two-cluster, horizontal mesh arrangement shown in Figure 7. 
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worked models and cluster material in the hot worked models. 
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Figure 38: Mesh of discrete model for large particles with a local volume fraction of 
4%, 
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Mesh of discrete model for small particles with a local volume fraction of 
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Mesh of discrete model for small particles with a local volume fraction of 


2%. 
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Figure 41: Mesh of discrete model for small particles with a local volume fraction of 
2%. Note that this mode! has circular shaped cluster regions. 
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The empty space surrounding the particles in Figure 37 is for visualization purposes only. 

The space in the models 1s filled in with matrix material in the cold worked models and 
~ 

cluster material properties in the hot worked models. The local volume percentage of 

particles in Figure 37 1s 25%. The particles were modeled as purely elastic material with 

Young’s modulus, E, equal to 60*10° DSI. 

The local volume fraction and size of reinforcement particles were also varied in 
the following manner. The reinforcement particle local volume fraction was reduced from 
25 % to 4% using the same particle size as shown in Figure 38; the particle size was 
reduced to one-quarter the size shown In Figure 37 and modeled with a local volume 
fraction of 4% and 2% (note Figures 39 to 41). 

B. STRESS AND STRAIN RESULTS OF THE DISCRETE MODELS 

The fringe plots of effective von Mises stress and effective plastic strain are shown 
in Appendix C. The trend of the results for all discrete models is similar. The discrete 
particles sustain the highest stresses and the lowest effective strains. Figures 42 and 43, 
for the 25% particle volume fraction and large particle size, apply to the low strain-rate, 
hot-work model. The stresses in the particles reach 875 ksi while the stresses in the 
matrix are as low as 15 ksi. The particles do not deform plastically while the effective 
plastic strains in the matrix are as high 3.60. The highest strains are found between the 
particle clusters, and above and below the particle clusters. Figures 44 and 45 show the 
results of the intermediate strain-rate, hot-work model with a particle size one-quarter that 
shown in Figures 42 and 43 and a volume loading of 4%. The results are qualitatively 


similar to those of Figures 42 and 43. 
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working at a low strain rate to a true strain of 1.20 with local particle 
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Figure 43: 


working at a low strain rate to a true strain of 1.20 with local particle 


volume fraction of 25%. 
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Figure 44: 


working at an intermediate strain rate to a true strain of 1.20 with local 


particle volume fraction of 4%. 
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Figure 45: 


working at an intermediate strain rate to a true strain of 1.20 with local 


particle volume fraction of 4%. 
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The horizontal displacement of the center of particles was measured for the 


different methods of TMP. These are shown in Figure 46 for the large particle size and 25 


™ 


volume pct., and in Figure 47 for the small particle size and 2 volume pct. 
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particle, 25% volume fraction, discrete model 






tone wa) ole’ ahs atwlais aatece tate’ 


ere 









law cate ote 
























anh path ty Se ah ag ho 














Snr eee 


chet A Rae a! A EsNehetet state "et Pat of dak ated A gate 






Displacement 








steht vhen 7 en weet tae dae atl cnet stp hate cite! el vha etd athe’ «atata't dhl wal rhe nscale 








Figure 46: Horizontal displacement of the center of particles for the large particle, 
25% volume fraction, discrete model. 
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igure 47: Horizontal displacement of the center of particles for the small particle, 
2% volume fraction, discrete model. 
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Figure 48 shows the height-to-width ratios, after the different methods of TMP, of 


the cluster regions shown in Figure 41. 


Height to width ratio after deformation for the two cluster discrete model with 
a local volume fraction of two percent 
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Figure 48: Height to width ratio after deformation of the small particle, 2% volume 


fraction, discrete model. 


Again, Figures 46 to 48 support the observation that the particles will move more 
readily using a hot vs. cold TMP. Figure 49 depicts a visual view of the cluster region 
deformation and particle movement for the small particle, 2% local volume fraction 
discrete model under different TMP methods. 

G. COMPARISON OF THE SMEARED VS. DISCRETE MODELS 

The smeared and discrete models do not respond in the same manner. First, the 

cluster zones have slightly different shapes in both models: the smeared model used a 


circular shaped cluster zone while the discrete model used a straight line approximation 


iscrete 


d 


10n 


he small particle, 2% local volume fract 


ion int 


* 


model after different methods of TMP to a true strain of 1.20. 


Visual view of cluster deformat 


— 


= 





= a SS 
— 


Cold worked 


— 
barnes 
—_ =" 


| 
my ae 


r 
— 
a 
= 
— 
= 
| 


~ —_— 
—— 
—_———; 
ad 
—— 
it 


ro 


gh strain rate 


h 


9 





Hot worked 





ih 


i 


= 


late strain rate 


in rate 


d 





n 


Visual view of cluster deformation in the small particle, 2% local volume 
fraction, discrete model with circular cluster regions after different methods 


of TMP to a true strain of 1.20. 
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Figure 49 


for a circle of the same size, as shown in Figure 37. One exception to this is the mesh 
shown in Figure 41. In the cold worked models the material properties of the particles 
were smeared into the cluster regions of the smeared model but ies the discrete model. 
On the other hand, the stiffness of the particles was neglected in the cluster zone of the hot 
work, smeared models. In the discrete models the stiff particles were modeled in both the 
hot and cold worked models. 

A qualitative comparison was made between the smeared and discrete models. 
When the horizontal displacement of the cluster center is compared as seen in Figures 33, 
46 and 47, it is noted that both models show the same general trend of greater 
displacement for the hot work models vs. the cold work models. The general trend of 
more displacement for the low strain rate vs. the high strain rate hot worked model is also 
seen. In Figure 47, it is seen that the hot work, low strain-rate model provides the largest 
amount of relative movement among the discrete particles. Again, this indicates that the 
best TMP to reach homogeneity, in the least strain is the hot work, low strain-rate method 
for both the smeared and discrete models. This finding 1s also supported by Figure 49. In 
Figures 27 and 48, which show the height-to-width ratios of the cluster regions after 
deformation, it is seen that the results from the discrete model are similar to those of the 
cluster models. The exception is that the results of the hot work model at high strain rates 
is not consistent. This case requires further investigation. In the smeared models it was 
assumed that larger deformation of the cluster region would result in more homogeneity. 
This assumption was supported from the discrete models as seen in Figure 49. This figure 


shows that for larger deformation of the cluster zone, the particles are more fully 


64 


dispersed. All of these comparisons give credibility to the work in Chapter ITI where the 
particle clusters are being modeled with smeared material properties. The smeared models 
are more efficient from a computational point of view and are considered a valid means of 
studying particle redistribution in a simplistic manner. 

An observation 1s made that with a higher local particle volume fraction in the 
discrete models, particle redistribution does not as readily occur as compared to a lower 
local particle volume fraction. This is illustrated visually in Figure 50. This leads to the 
observation that, if the particle clusters in the consolidation process are tightly packed, 
particle homogeneity is harder to achieve. If the particle clusters are loosely packed, 1.e., 
better mixing in the consolidation/fabrication technique, then particle homogeneity will 
more easily be reached. If the cluster zone has many particles, the particles will support 
the most load in the cluster zone and dominate the behavior of the cluster deformation. 


As a result the particle redistribution appears to be less likely. On the contrary, the matrix 


dominates the cluster zone behavior with a low local particle volume fraction. 


Ing atvan 


fle 2Os 


‘ete models after hot work 


ISCI 


in the d 
late strain rate to a true strain O 


‘mation 


“med 


intel 


Visual view of cluster defo1 



















ATW 
Te fl Y i 
Hi ih 
Hi : Hh 
he inh 
wi iy 
Hi tT 
h 
! Hh ‘ 
i | 
rd | 
= Ne 
r N 
t Mi HI | 
AYR ih 
ee HTH 
© tt 
~~ = TN] 
O oO O nH 
= © = , 
oO OQ i 
o ; = t iH 
— ~ © I in 
2 D 3 i 
te a O HTT 
a n ag 
5 | be ie HH 
— { ce | 
< i | 2 = k 
ss mu T| = Oo TOME OSU 
y tg int eae 
= He en “S aR 
= = i | Tn SO Kf 
= Ht wi mt = CMT 
© H = Nea ial = 7 na Note yuna 
e Ee Hi O H = a Mt af 
4 em r > : 
3 ee E itt 
2 sa 
a c i) ie S rH 
aM) : H 3. e i: 
HIATT 1 fl ~ ° ; 
iN t a : HS s ui 
Me Ta [ H oO n rt o im AL 
hit nt AIT — v il 
A : 3 2 it 
We HN C ca = UL 
AAV LD E tC th 
ni WH ae wo 1 | 
a HE Oo. Oo. il nT 
Pe ep » = ait 
be! ae a a = i i 
_ ‘Ss F xl 


Oo 
o 


felez 0: 


in the discrete models after hot workin 


1i0n 
te to a true Strain oO 


in Ta 


66 


f cluster deformat 


diate stra 


Visual view oO 
at an interme 


° 
* 


Small particle, 2% local volume fraction discrete model 


Figure 50 


V. SUMMARY 


AS. CONCLUSIONS 

The following conclusions are drawn from this work: 

l. The current modeling methods, both smeared and discrete, used to model 
MMC deformation is an inexpensive and informative means of evaluating the matrix and 
cluster material behavior during the entire deformation process. 

a After cold working all the smeared models exhibited higher stresses and 
lower effective plastic strains within the cluster regions. 

Se After hot and cold working all the discrete models exhibited higher stresses 
and lower effective plastic strains in and around the particles. 

4. After hot working at the low strain rates, all the smeared models exhibit 
lower stresses and higher effective plastic strains within the cluster region as compared to 
the matrix region. 

S, After hot working at the intermediate strain rates all the smeared models 
exhibit higher stresses and higher effective plastic strains within the cluster region as 
compared to the matrix region. These higher stresses in the cluster region only become 
apparent near the end of the deformation process and is a result of the fact that the true 
strain rate is higher near the end of the hot working process compared to the beginning of 
the hot working process. 

6. After hot working at the high strain rates all the smeared models exhibit no 


distinction between the cluster regions and the matrix regions. 


67 


dl. The results, using the linear law compared to the quadratic law, in the 

cluster regions of the smeared models after hot working at intermediate and low strain 
~ 

rates, are as follows. The difference in effective plastic strains between the linear and 

quadratic laws is very small for both strain rates. The linear law results in higher stresses 

within the cluster regions for the intermediate strain rate model and lower stresses for the 

low strain rate model, when compared to the quadratic law. 

8. Internal stresses achieved in the cold worked models are an order of 
magnitude higher than those achieved during the hot worked models. 

9, Internal stresses achieved in the hot worked, high strain rate models are 
twice the value as those found in the hot worked, low strain rate models. 

10. Particle homogeneity within the matnx material 1s better achieved using a 
hot working process vs. a cold working process. This is supported by the deformability of 
the clusters in the smeared models as measured in terms of height to width ratios of the 
clusters after straining, the vertical compression of the cluster material compared to the 
matrix material after straining and the horizontal displacement of the cluster centers during 
the straining process. The relative particle movement tn the discrete models also support 
this finding. 

11. ‘Particulate homogeneity within the matrix material is better achieved using 


a lower strain rate hot working process. This is supported by the same measures 


described in the previous conclusion. 


68 


B. RECOMMENDATIONS 

The following recommendations are made for subsequent studies: 

ie The current model does not allow breaking down of a single cluster into 
smaller clusters and mixing within the surrounding matrix material. A further experimental 
study 1s needed to gain more insight of the actual physical occurrence within MMC’s 
during TMP. The subsequent study should improve the modeling technique to more 
closely describe the actual phenomena. For this description the Eulerian-Lagrangian 
technique of the FEM might be more suitable. 

Ze The current FEM package used does not allow for combining constitutive 
equations of a material region for different mechanisms of creep. During the deformation 
process, the mechanisms of creep can change based on current stresses or strain rates. It 
is suggested that in the further development of this study, additive constitutive laws pod 
to be incorporated into the FEM program. 

By The discrete model needs to be further developed. Focus should be on the 
effects of initial particle distribution within the cluster regions including size, arrangement 
and local volume fraction, and the size or boundary of the cluster zones. Also the 


discrepancy between the smeared model and discrete model for the hot worked high strain 


rate case needs further investigation. 
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APPENDIX A: FRINGE PLOTS OF EFFECTIVE VON MISES STRESS AND 
EFFECTIVE PLASTIC STRAIN FOR ALL COLD WORKED SMEARED 
MODELS. 
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Figure A2 


arrangement in the cold worked model at a true strain of 1.20. 
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APPENDIX B: FRINGE PLOTS OF EFFECTIVE VON MISES STRESS AND 
EFFECTIVE PLASTIC STRAIN FOR ALL HOT WORKED SMEARED 
MODELS. 


81 





LEO t/t ‘6 'S ‘6h 


So Zee 


time = 3.49998E-G2 
fringes of eff. stress (v-m) 


min= 6.333E+@3 in element 822 
max= 1.@51E+84 in element 572 





Or 


fringe levels 
7.7A8E+83 5 


7.956E+83 ™ 
B8.211F+03 > 


B.467E+83 ° 


~—— 9. 722F+83 


8,2393E+83 > 


9.409E+83 : 


9. 744E+89 > 


1, GUBE+B4 





Fringe plot of Effective Von-Mises stress for the two cluster horizontal 


Figure B1 


arrangement in the hot worked high strain rate model at a true strain of 


1.20. 


82 








x 

— me 
s « 
ab 

S 2 
T" | | 
3 3 
yo o 


Figure B2: 


LG/BK/b “PP'9 ‘6t 


©) 
33 
me “h ck a 
x 3 > oe, 
| a 4 a 
mM 
N = ca a 
. 2. fo tt 
“iN a i 
og TU 
NN Oo WwW N 
m TT + 
+ + Dd 
Q ad m W 
8 & +h Ww 
+h Ww 
me oe . @M 
Sees Mm 
ce I 
0 — © 
— ome »N 
m 0 
334 Ca 
o © ae 
es QO 
ce ce 
0 
Mme ce 
“NW os 
NN wy 
oY mae 
Ba | 
§ 
WES 
AT ETS 
SOU 
+4 
; % 
= — et *] ~ rs be m > 
w ay 3 i) w Ww) w w CO 
© a “I =) 07 n a vom 
2D Ww eee nN “J Re a Qo 
mM ™ Mt ™m Mm rl ™ ™m 
+ + + + i = + + + = 
oS So a 3 © o © S 
ao & 8 OE & ao g a a 
=r mM ex 
— 
a =| 
Oo 





Fringe plot of Effective plastic strain for the two cluster horizontal 
arrangement in the hot worked high strain rate model at a true strain of 


20: 


83 








96/6 /EN'ZE-PL El 


©) 


Ti WHE H Hl nh WH 5 3 Se WN 
N << walk i i tH ait s 3 2 5 _ 
we a 
‘ i 2 e Mo U U 
y ni ee 
tt 
HALA s eu . ww (WW) 
ae ae a = 
i TL q 
t UHI (PN w Ww ice 
Be -—2 fee e Ww 
Sa mm 
“” | 
cv ae) 
— os le 
vc % 
33 ” 
o “” 
a = 
ce ck 7 
< 
pe AJ t 
ng 3 
Ww N ~ 
«A 
~h 
| 
> - ” ” Y ui a ul oes 
tw 7 ;: . 
2 q a 0 wo a x ur w ca 
Q © Ov Ww w a +b N Qo 
w) Ww w w w a w w w } 
o 
a 





Figure B3: _—‘ Fringe plot of Effective Von-Mises stress for the two cluster horizontal 
arrangement in the hot worked intermediate strain rate model at a true 


strain of 1.20. 


84 





46/6 /b LE Lect 


(9: 









eH 
7 Ol me h ce a 
re bts} it 3 4s) x J > | = 
La (ies cea a 
srounet tH ti ii a HDT Hi N ca “UD 
Au SL Hi ah saseite a ey) Z 4s fee 
eT pat fH ea vl @ a 
+t)! A ie i HH Hy iO N“N Go Ww {) 
f Lit rH a Sa 
Hae ith eth + + > 
a) Q a m W 
Q a h 
h 
oe mat » W 
= | m™ 
oT I 
0 © —- @ 
bt tee fv 
9 i) 
a 3 ct 
2 0 
el & 0 
oe ct 
V 
Oe + 
“NO 74 
NN wu 
Wo 
3 





= 


P a co 
& . c , e 
lee Te ciate ace th ae ce ated “hy 
- ee 


i Ty ats | TDL 

{ 

MEAS (ATLL Ree Ture 
fakes) Hs 


+ vba M 
7H rt | 
Lt 


J] + 
Aid | 


66+362+r ‘Tt 
BB+3+Tr “Tt 
@6+366¢ “tT 
88+3EG¢r ‘Tt 





Oo+IZ66 T 

@0+348E "Tt 
§ @B+3TGE ‘t 
; OB+392E°T 





B8+342E°T 
s[anay aBuru3 





Figure B4: ‘Fringe plot of Effective plastic strain for the two cluster horizontal 
arrangement in the hot worked intermediate strain rate model at a 


true strain of 1.20. 
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Figure B6: _—‘ Fringe plot of Effective plastic strain for the two cluster horizontal 
arrangement in the hot worked low strain rate model at a true strain of 
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Figure B7: _‘ Fringe plot of Effective Von-Mises stress for the two cluster vertical 
arrangement in the hot worked high strain rate model at a true strum o! 
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Figure B9: _—‘ Fringe plot of Effective Von-Mises stress for the two cluster vertical 
arrangement in the hot worked intermediate strain rate model at a true 


strain of 1.20. 
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Figure B10: Fringe plot of Effective plastic strain for the two cluster vertical 
arrangement in the hot worked intermediate strain rate model at a true 


strain of 1.20. 
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Figure B11: Fringe plot of Effective Von-Mises stress for the two cluster vertical 
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Figure B12: Fringe plot of Effective plastic strain for the two cluster vertical 
arrangement in the hot worked low strain rate model at a true strain of 
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Figure B13: Fringe plot of Effective Von-Mises stress for the square cluster 
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Figure B14: Fringe plot of Effective plastic strain for the square cluster arrangement In 
the hot worked high strain rate model at a true strain of 1.20. 
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Fringe plot of Effective Von-Mises stress for the square cluster 
arrangement in the hot worked intermediate strain rate model at a true 


strain of 1.20. 
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Figure B16: Fringe plot of Effective plastic strain for the square cluster arrangement in 
the hot worked intermediate strain rate model at a true strain of 1.20. 
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Figure B18: Fringe plot of Effective plastic strain for the square cluster arrangement in 
the hot worked low strain rate model at a true strain of 1.20. 
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Figure B19: Fringe plot of Effective Von-Mises stress for the square cluster 
arrangement in the hot worked intermediate strain rate model using the 


linear law at a true strain of 1.20. 
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Figure B20: Fringe plot of Effective plastic strain for the square cluster arrangement in 
the hot worked intermediate strain rate model using the linear law al a truc 


strain of 1.20. 
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Fringe plot of Effective Von-Mises stress for the square cluster 


Figure B21: 


arrangement in the hot worked low strain rate model using the linear law 


at a true strain of 1.20. 
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Figure B22: Fringe plot of Effective plastic strain for the square cluster arrangement in 
the hot worked low strain rate model using the linear law at a true strain of 
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Figure B23: Fringe plot of Effective Von-Mises stress for the diamond cluster 
arrangement in the hot worked high strain rate mode! at a true strain of 
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Figure B26: Fringe plot of Effective plastic strain for the diamond cluster arrangement 
in the hot worked intermediate strain rate model at a true strain of 1.13. 
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Figure B28: Fringe plot of Effective plastic strain for the diamond cluster arrangement 
in the hot worked low strain rate model at a true strain of 1.13. 
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APPENDIX C: FRINGE PLOTS OF EFFECTIVE VON MISES STRESS AND 
EFFECTIVE PLASTIC STRAIN FOR ALL DISCRETE MODELS 
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Fringe plot of Effective Von-Mises stress for the discrete cold worked 
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Figure C2: _ Fringe plot of Effective plastic strain for the discrete cold worked 
model, with a local volume fraction of 25%, at a true strain Otel: 
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Fringe plot of Effective Von-Mises stress for the discrete hot worked high 
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Figure C4: Fringe plot of Effective plastic strain for the discrete hot worked high 
strain rate model, with a local volume fraction of 25%, at a true strain of 
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Figure C15: Fringe plot of Effective Von-Mises stress for the discrete hot worked 
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Figure C17: Fringe plot of Effective Von-Mises stress for the discrete hot worked low 
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Figure C19: Fringe plot of Effective Von-Mises stress for the discrete hot worked 


intermediate strain rate model, with a local volume fraction of 4%, ata 


true strain of 1.20. 
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